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1. Introduction 

The origin of the nonthermal continuum emission in Blazars is a challenging 
problem in the area of active galactic nuclei. The X-ray emission is usually 
interpreted as synchrotron radiation of relativistic electrons accelerated in 
a jet which itself moves at relativistic speed towards the observer, while the 
TeV gamma-rays are though to arise via synchrotron self-comptonisation 
(e.g. Mastichiadis & Kirk 1997), via comptonisation of external photons 
(e.g. Dermer, Schlickeiser, Mastichiadis 1992) or via the decay of relativis- 
tic protons (e.g. Mannheim 1998). 

In this contribution, we present a self-consistent model, where electrons ex- 
perience first-order Fermi acceleration at a shock front in a relativistic jet 
and radiate synchrotron emission in a post-shock region which contains a 
homogeneous magnetic field. The full time, space and momentum depen- 
dence of the electron distribution function is used for a calculation of the 
nonthermal synchrotron spectra. By varying the rate at which electrons 
enter the acceleration process, we calculate the time-dependence of the 
spectral index for the X-ray emission predicted in this model. The results 
show that the synchrotron spectral index displays a characteristic looplike 
behaviour with intensity (as has been observed in several blazars), where 
the orientation of the loop depends on whether the acceleration time scale 
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is comparable to the synchrotron cooling time scale or not. We show that 
our model provides a good fit to the observed evolution of the spectral 
index of Mkn 421 during a flare in 1994. 

2. Acceleration model 

We assume diffusive particle acceleration at a shock front moving at a con- 
stant speed Us along a cylindrically symmetric jet. Following Ball &: Kirk 
(1992) we divide the jet into an acceleration region close to the shock and a 
cooling region downstream (see figure). Electrons in the acceleration region 
repeatedly encounter a plane shock. We assume that high energy particles 
pass unaffected through the shock and scatter elastically of magnetic irreg- 
ularities (Alfvcn waves) in the fluid on either side of the sliock. In the rest 
frame of the shock the speed of the upstream fluid exceeds the downstream 
fluid speed so that particles gain energy by the first-order Fermi process 
(Drury 1983). Electrons which escape from the acceleration region into the 
downstream region cease to diffuse and are assumed to remain frozen into 
the fluid flow, suffering only synchrotron losses. 

x = 




upstream acceleration zone downstream radiating region x 

Consider electrons which are injected and enter the acceleration mecha- 
nism at the constant rate Q with Lorentz factor 79. Although it is possible 
to write down the distribution function of such particles as a function of 
momentum, space and time (e.g. Toptygin 1980, Drury 1991), we use a 
spatially averaged model (Axford 1981, Bogdan Sz Volk 1983), in which an 
electron is presumed to undergo continuous acceleration while in the vicin- 
ity of the shock, such that its Lorentz factor 7 increases at a rate 7/iacc) 
where = D/tc is the acceleration rate and D is the average fractional 
momentum gain per cycle (for crossing and recrossing tlic shock), and tc is 
the average time taken to perfom such a cycle (sec Bell 1978, Drury 1983, 
Kirk, Melrose, Priest 1994). The acceleration process is coupled with the 
fact that particles in the downstream region have a chance of not return- 
ing to the shock. These particles escape at a rate which is just the 
escape probability per cycle Pesc divided by the average cycle time tc- As- 
suming pitch angle scattering is sufficiently rapid to keep the distribution 
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function almost isotropic, we can easily include the electron synchrotron 
loss term (7)5 = —/3s 7^ in a differential equation governing the number 
A'^(7)d7 of particles in the acceleration zone around the shock with Lorentz 
factor between 7 and 7 + d7 (e.g. Kirk, Melrose, Priest 1994, Mastichiadis 
& Kirk 1997): 



dN d_ 



- Ps 7') ^1 + 7^ = Q^ii - 70) (1) 

where Ps = 4: / {6 fiQ nie c) with (Tt = 6.65 • 10~^^m^ the Thomson 
cross section and B the magnetic field in Tesla. There obviously exists an 
important Lorentz factor in this problem: at the point 7max = l/(/?s^acc) 
acceleration is exactly balanced by synchrotron losses. It is to be expected 
that the solution always vanishes for 7 > 7max and also for Lorentz factors 
to which there has not been enough time to accelerate particles. 
Prom the synchroton losses we get a typical cooling time scale given by 
icooi = 1/(7 /3s). As 

7 — ^ 7max) ^cool tcnds to tacc- 

In the model developed by Ball &: Kirk (1992), accelerated particles es- 
cape into the downstream plasma, where they radiate. We can formulate 
the kinetic equation obeyed by the density of particles in the radiation 
zone most compactly using a coordinate system at rest in the radiating 
plasma. Note that in this case the radiation must be Doppler-boosted to 
the observer's frame: quantities such as the intensity or the flux of the 
synchrotron radiation, which are computed in the rest frame of the source 
(plasma rest system), require a Lorentz transformation into the observer's 
frame, determined by the Doppler factor 6 = l/[r(l — P cos6)], where 
P = (1 — M^/c^)"^/^, P = U2/C and 9 is the angle between the direction of 
motion and the line of sight to the observer. 

In the plasma rest system the shock front provides a moving source of elec- 
trons, which subsequently suffer energy losses, but are assumed not to be 
transported in space. The kinetic equation governing the differential density 
dn{x,j,t) of particles in the range dx, dj then is 

where Xs{t) is the position of the shock front at time t. Equations (1) and 
(2) have been solved analytically (see Kirk, Rieger, Mastichiadis 1998). 
To obtain the synchrotron emissivity as a function of position, time and 
frequency, it remains to convolve the density n with the synchrotron Green's 
function P(i/, 7) (e.g. Melrose 1980). Then at a point x = X (> Ust) on the 
symmetry axis of the source at time t, the specific intensity of radiation in 
the X direction depends on the retarded time t = t — X/ c and is given by 

I{u,t) = j d7P(z/,7) j dxn(x,7,t + x/c) (3) 
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3. Spectral variability 

The model described above provides a remarkably good fit to the radio 
to X-ray spectra of several blazars. In a recent publication (Kirk, Rieger, 
Mastichiadis 1998) we apply our model to the nonthermal X-ray emission 
of Mkn 501, where we get a very good representation of the observed syn- 
chrotron spectrum. Since we solve Eq. (1) and Eq. (2) for the full time, 
space and momentum dependence of the electron distribution function, we 
arc able to apply our model to the observed variability or flaring behaviour 
in blazars. Variability or flaring can arise for several reasons, e.g. when the 
shock front overruns a region where the local plasma density is enhanced. 
In this case the number of particles entering the acceleration process might 
increase for a time tj. A simple representation of a flare, for example, is 
obtained by setting Q{t) = Qq foi t < and t > ti and Q{t) = Qi for 
< t < ti- In Fig. (1) we plot the corresponding synchrotron emission 
given in the rest frame of the downstream plasma for the case Qi = 1.8 Qo- 
Well below the maximum emitting frequency, where the acceleration time 
scale tacc is much shorter than the synchrotron cooling time scale tcooi) the 
spectrum shows the typical soft lag behaviour observed in several sources 
(e.g. PKS 2155-304: Sembay et al. 1993, H0323-022: Kohmura et al. 1994, 
Mkn 421: Takahashi et al. 1996): the spectrum becomes softer in declining 
phases and harder during phases of rising flux, which indicates that vari- 
ations in the hard X-rays always lead those in the soft X-rays. An inter- 
pretation of this behaviour has been given in terms of synchrotron cooling 
(e.g. Tashiro et al. 1995). In this case, information about changes in the 
injection propagates from high energy, where the cooling process is faster, 
to low energies. 

In Fig. (1) we also give a qualitative fit for the evolution of the photon in- 
dex a for Mkn 421 during an X-ray flare, observed with the X-ray satellite 
ASCA in May 1994 (note that the data should be treated with caution, 
since they were derived for the physically unlikely case that the absorbing 
column Nh was allowed to vary, see Takahashi et al. 1996). For the plot 
we use both the constraints given in Takahashi et al. (1996) and the pa- 
rameters obtained by Mastichiadis & Kirk (1997) in fitting the observed 
multiwavelength spectrum of Mkn 421 in 1994. The observed maximum 
frequency is roughly fmax ~ 1.4- lO^^Hz, so that, applying a Doppler factor 
of 15, the maximum electron Lorentz factor is 10®. The results are plotted 
in Fig. (1) for a shock velocity Ug = 0.26 c and for frequencies of obser- 
vation well below the the maximum frequency. Closer to the maximum 
frequency, where the acceleration and cooling time scale become compara- 
ble, the picture changes. In Fig. (2) we show the predicted results for the 
same parameters as used in Fig. (1), for a frequency of observation close to 
the maximum frequency. Here the loop in the photon index a is traced in 
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Figure 1. Evolution of the photon index a and intensity (arbitrary units) in a flare at 
low frequencies. The loop in the photon index a is followed in the clockwise direction. 
The Mkn 421 data points in this picture are taken from Takahashi et al. 1996 (see text 
for details). 

the anti-clockwise direction. This behaviour is a typical signature of the un- 
derlying acceleration mechanism. As we have argued recently (Kirk, Rieger, 
Mastichiadis 1998), this feature represents the effect that now information 
about the occurrence of a flare propagates from low to high energy as par- 
ticles are gradually accelerated. 

Although in our model we neglect factors such as, for example, inhomogen- 
ities in the source geometry, the result that the loop is traced in a clockwise 
direction if controlled by synchrotron cooling and followed in the opposite, 
i.e. anti-clockwise direction if dominated by the acceleration time scale is 
robust. 
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Figure 2. Evolution of the photon index a and intensity (arbitrary units) in the same 
flare as plotted in Fig.(l) if the system is observed near the maximum frequency z/max. 
The loop in the photon index a is now followed in the anti-clockwise direction. 
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